The stochastic dynamics of micron and nanoscale elastic cantilevers in fluid: 

fluctuations from dissipation 
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The stochastic dynamics of micron and nanoscale cantilevers immersed in a viscous fluid are 
quantified. Analytical results are presented for long slender cantilevers driven by Brownian noise. 
The spectral density of the noise force is not assumed to be white and the frequency dependence 
is determined from the fluctuation-dissipation theorem. The analytical results are shown to be 
useful for the micron scale cantilevers that are commonly used in atomic force microscopy. A 
general thermodynamic approach is developed that is valid for cantilevers of arbitrary geometry as 
well as for arrays of multiple cantilevers whose stochastic motion is coupled through the fluid. It 
is shown that the fluctuation-dissipation theorem permits the calculation of stochastic quantities 
via straightforward deterministic methods. The thermodynamic approach is used with deterministic 
finite element numerical simulations to quantify the autocorrelation and noise spectrum of cantilever 
fluctuations for a single micron scale cantilever and the cross-correlations and noise spectra of 
fluctuations for an array of two experimentally motivated nanoscale cantilevers as a function of 
cantilever separation. The results are used to quantify the noise reduction possible using correlated 
measurements with two closely spaced nanoscale cantilevers. 



I. INTRODUCTION 



The dynamics of micron and nanoscale cantilevers are 
important to a wide variety of technologies. For example, 
the invention of the atomic force microscope (AFM) yij , 
which relies upon the dynamics of a cantilever a few hun- 
dred microns in length, has revolutionized surface science 
paving the way for direct measurements of intermolecu- 
lar forces and topographical mappings with atomic pre- 
cision for a broad array of materials including semi- 
conductorSjDolymers, carbon nanotubes, and biological 
cells @,0,yriSlH0 ( see IE 13 for current reviews). In 
conventional dynamic atomic force microscopy the can- 
tilever is used to measure the force interactions between 
the cantilever tip and sample. Cantilevers smaller than 
conventional AFM have also been used to unfold single 
protein molecules with improved force and time resolu- 
tion (Tjj. I* nas a l so been proposed to exploit the inher- 
ent thermal motion of small cantilevers to make dynamic 
measurements of single molecules [l2j . A passive un- 
driven cantilever placed in a viscous fluid will exhibit 
stochastic oscillations caused by the thermal bombard- 
ment of fluid molecules by Brownian motion. In fact, 
measuring the thermal spectra of the cantilever is a com- 
monly used AFM calibration technique [l3L ITU . 

In all of these applications the ultimate force sensitiv- 
ity of a particular measurement is limited by the inherent 
thermal noise of the experimental system. There are at 
least two ways to improve upon this limitation: (i) make 
the cantilevers smaller [Tfjl ITU , (ii) use correlated mea- 
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surements of multiple cantilevers 0, 0] . 

Uniformly decreasing the dimensions of a cantilever 
results in the favorable combination of decreasing the 
cantilever's equivalent spring constant while increasing 
its resonant frequency yielding improved force sensitivity 
and time resolution. By measuring the cross-correlations 
between two cantilevers in fluid the independent fluc- 
tuations of the two do not contribute, leaving only the 
smaller correlated fluctuations due to coupling through 
the fluid. This type of approach has been used to mea- 
sure femtonewton forces on millisecond time scales be- 
tween two micron scale beads placed in separate optical 
traps (an improvement of a hundredfold from prior mea- 
surements) [l7j . The ability to significantly increase the 
force resolution by making correlated measurements has 
yet to be exploited for micron and nanoscale cantilevers. 
The magnitude of the fluid coupled noise will depend 
upon the spacing and exact geometries of the cantilevers. 
Combining (i) and (ii) and measuring the correlations of 
multiple nanoscale cantilevers offer the potential for ex- 
perimental measurements with unprecedented force and 
time resolution (TUTIITH. 

As experimental measurement continues to push to- 
ward the stochastic limit it is important that we build 
a physical understanding of the stochastic dynamics of 
micron and nanoscale cantilevers for the precise condi- 
tions of experiment including complex cantilever geome- 
tries IHIll, the effects of nearby walls [H HJ lH E3 , 
and the fluid-coupled d yna mics of multiple cantilevers in 
an array configuration 12]. 

Although micron and nanoscale cantilevers exhibit 
stochastic motion due to the thermal motion of matter 
the elastic structures are still large compared to indi- 
vidual fluid molecules and the equations of continuum 
mechanics remain valid. In what follows we are con- 
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cerned with situations where the Knudsen number (the 
ratio of the mean free path of the fluid particles to the 
width of a cantilever) remains sufficiently small so that 
this statement remains true. This means that the fluid 
is described by the usual Navier-Stokes equations with 
no-slip and stress continuity boundary conditions at the 
solid surfaces. 

The Navier-Stokes equations governing the motion of 
an incompressible fluid, and written in nondimensional 
form, are 



du -> _, 

at 



-Vp + V 2 u, 



0, 



(1) 
(2) 



where u is the fluid velocity, p is the pressure, and t is 
time. There are two inertial terms on the left hand side 
of Eq. Q multiplied by the nondimensional parameters 
Ruj and R u . The Strouhal number R u — L 2 jvT plays the 
role of a frequency based Reynolds number expressing the 
ratio between local inertia forces and viscous forces where 
L and T are characteristic length and time scales, respec- 
tively. The velocity based Reynolds number R u = UL/v 
expresses the ratio between convective inertial forces and 
viscous forces. 

Micron and nanoscale cantilevers are characterized by 
high oscillation frequencies and small oscillation ampli- 
tudes. In this case R u C 1 so that the nonlinear convec- 
tive inertial term tt«Vtt is negligible and the equations 
become linear. However, the Strouhal number, 
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(3) 



(using the half- width w/2 as the appropriate length scale) 
is often not negligible. As a result, the local inertia term 
must be kept in Eq. Q making the analysis more diffi- 
cult. In addition, experimentally motivated cantilevers 
are often of complex geometry, near surfaces, or in an ar- 
ray configuration with multiple cantilevers in close prox- 
imity. These difficulties have led to the development of a 
thermodynamic approach to calculate the stochastic dy- 
namics of micron and submicron scale cantilevers for the 
precise conditions of experiment (discussed below) |12| . 
In the following it is assumed that R u is negligible, and 
we leave off the subscript denoting R w by R. 

The fluid equations are coupled with the governing 
equations of elasticity, 



d 2 w 



V.cr, 



(4) 



where a is the stress tensor, p c is the cantilever den- 
sity, and im(x, y, z, t) is the cantilever deflection. Equa- 
tions C3),@, and Q represent the governing determinis- 
tic continuum equations. If one considers lengths scales 
on the order of a few atomic lengths these equations 
would include stochastic terms 23]. This represents 
a difficult fluid-solid interaction problem governing the 
stochastic cantilever dynamics. A brute force molecular 



dynamics approach which would resolve the stochastic 
motion of all of the fluid and solid molecules is computa- 
tionally prohibitive. However, the system is always near 
equilibrium permitting a much more accessible thermo- 
dynamic based solution strategy. 

The paper is organized as follows. A very general 
thermodynamic approach, valid for complex geometries 
as well as for arrays of closely spaced cantilevers, is 
discussed in Section |H] Analytical theory based upon 
the thermodynamic approach, valid for long slender can- 
tilevers, is discussed in Section IIIII In Section IIVI a mi- 
cron scale cantilever, similar to an atomic force micro- 
scope, is explored and in Section [V] the stochastic dy- 
namics of an array of experimentally motivated nanoscale 
cantilevers are quantified using the thermodynamic ap- 
proach with deterministic finite element numerical sim- 
ulations to generate results for the precise conditions of 
experiment. A general discussion of the approach can be 
found in Ref. 
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FIG. 1: Schematic illustrating an array of two opposing can- 
tilever beams separated by a distance s. In the thermody- 
namic approach of Section [21 a step force F\ is removed from 
the cantilever on the left and the deterministic cantilever de- 
flections Xi(t) and X2(t) are calculated. 



II. THERMODYNAMIC APPROACH 

The thermodynamic approach discussed here is based 
upon the fluctuation-dissipation theorem which states 
that for equilibrium systems the manner in which the 
system returns from a linear macroscopic perturbation 
is related to time correlations of equilibrium microscopic 
fluctuations Underlying this statement is the 

fact that equilibrium fluctuations and the dissipation of 
the system responding to a macroscopic perturbation are 
governed by the same physics. For the case of miniature 
elastic objects in a fluid the dissipation is dominated by 
the fluid viscosity and the fluctuations by the motion due 
to the bombardment by the fluid molecules. In what fol- 
lows, we assume that all of the dissipation comes from 
the fluid and that elastic dissipation in the cantilever is 
negligible. This assumption, however, is not required and 
other sources of dissipation could be included if desired. 
The essence of our approach is to calculate the dissipation 
deterministically and to use the fluctuation-dissipation 
theorem to determine the cantilever's stochastic dynam- 
ics. The approach is exact and the only assumptions that 
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have been made are classical dynamics and linear pertur- 
bations from equilibrium. The deterministic calculation 
of the dissipation can come from analytical theory, sim- 
plified models, or from detailed numerical simulations. 
The major benefit of this approach is that the determinis- 
tic calculations are straightforward, not computationally 
prohibitive, and methods of calculation are sophisticated 
and readily available. 

We will introduce the use of the thermodynamic ap- 
proach for the case of two opposing cantilevers as shown 
in Fig. 2] Consider one dynamical variable to be the 
stochastic displacement of the cantilever on the left X\ (t) 
where Xx(t) is a function of the microscopic phase space 
variables consisting of 37V coordinates, r 3N , and conju- 
gate momenta, p 3N , of the system and N is the number 
of particles in the system. The statistical treatment that 
follows determines the ensemble average of the cantilever 
deflections x%(t) over all experimental possibilities. For 
the undisturbed equilibrium system the ensemble aver- 
age, and correlations, are denoted as (). 

We now take the system to a prescribed macroscopic 
excursion from equilibrium and observe how the system 
returns to equilibrium. The macroscopic cantilever de- 
flection as it returns to equilibrium from a prepared ini- 
tial condition will be denoted as X\ (t) . For clarity of pre- 
sentation, lower case variables are reserved for stochastic 
quantities and upper case variables are for determinis- 
tic quantities. The connection between the stochastic 
and deterministic quantities is particularly straightfor- 
ward when the excursion from equilibrium is achieved 
through the application of a step force some time in the 
distant past that is removed at time t = 0, i.e. the force 
given by, 



/(*) 



Fx for t < 
for t > 0. 



(5) 



The applied force should be conjugate to the variable for 
which the fluctuations are to be calculated. For the can- 
tilever system under consideration here, to determine the 
fluctuations of the tip displacements, the force is applied 
to the tip of the left cantilever. The step force couples 
with the cantilever deflection Xx(t) and the full Hamilto- 
nian of the system H is given by, 



H = 



H + AH 
Ho 



for t < 
for t > 



(6) 



where AH = —FxXx(t) — — iqxi(0) and H is the un- 
perturbed Hamiltonian. For small perturbations, Fx, 
the quantity AH will also be small, and the results are 
greatly simplified. The nonequilibrium ensemble average 
Xx (t) is given by, 



Xx(t) 



J dr 3N dp 3N Xi ^ e -0(Ho+AH) 
J d r 3Nd p 3N e -f3(H +AH) ' 



(7) 



where j3 = (kgT) , kg is Boltzmann's constant, and 
T is the absolute temperature. In Eq. 0) the notation 



x\(t) is used to represent the value of xx evaluated at 
the phase space coordinates that have evolved from the 
values r and p at time t = |2f| . In the limit of linear 
perturbations, AH <C 1, this simplifies to, 

Xx(t) = ( Xl )+f3(AH) (xx)-(3(xxAH)+0(AH% (8) 

where an equilibrium ensemble average is given by, 

Jdr m dp m x 1 (t)e-> 3Ho 



(xx(t)} 



Jdr^dp^e-PH" 



(9) 



If we now assume that x\ and Xx have the equilibrium 
average (xi) subtracted and recall that AH = — i<\xi(0), 
we have our desired result, 



Xx(t)=/3Fx(xx(t)xx(0)) 
which can be rearranged to yield, 

(xx{0)xx{t)) = k B T 



Xx{t) 



(10) 



(11) 



The analysis is similar for the cross correlations of the 
deflections for two cantilevers in an array, 



(xx(0)x 2 (t)) = ^X 2 (t), 



(12) 



where X 2 is the displacement of tip 2 arising from the 
step force Fx applied to tip 1. 

It is interesting to highlight that the cantilever mo- 
tions are indeed correlated as indicated by Eq. 1)12(1 . The 
correlated fluctuations appear contrary to the naive idea 
that random molecular impacts upon the individual can- 
tilevers should only lead to uncorrelated motion. Ad- 
ditionally, it should be emphasized that the correlated 
motion of two (or more) cantilevers is no more difficult 
a calculation than determining the autocorrelation of a 
single cantilever. 

On the left hand side of Eqs. 1(11(1 and ((12(1 are stochas- 
tic quantities and on the right hand side are deterministic 
quantities. This permits a significant reduction in effort 
in the calculations of the auto- and cross-correlations of 
the cantilever deflections. The macroscopic deflections 
Xx(t) and X 2 (t) can be found from deterministic the- 
ory, model equations, or numerical simulations and the 
stochastic quantities are found simply through applica- 
tion of Eqs. ((11(1 and ((12(1 . The deterministic problem 
then reduces to the fluid-solid interaction problem given 
by the removal of a step force on a cantilever. For the 
step force Eq. 10), for t < 0, Xx(t) has a finite deflec- 
tion. After the removal of the step force the cantilever 
returns to its equilibrium position of Xx — 0. The precise 
manner in which the cantilever returns to equilibrium is 
dominated by the dissipation in the fluid. An adjacent 
cantilever will exhibit time dependent deflections given 
by X 2 {t). In this case both the initial and final equilib- 
rium deflections are X 2 = 0. The motion of the second 
cantilever is a result of the fluid motion caused by the 
first cantilever as it returns to equilibrium. 
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The spectral properties of the auto- and cross- 
correlations are determined by taking the cosine Fourier 
transform of Eqs. (f 1 1|) and Ijl2|l with appropriate factors 
(see Eq. $lbl ). This yields the noise spectra, Gn(w) and 
Gi 2 (w), given by, 

C? u (u) = 4 / (xi(0)a;i(t))cos(^)dt, (13) 
Jo 

/>OC 

G 12 (uj) =4 (xx(0)x 2 (t)) cos(ujt)dt, (14) 



where u> is the angular frequency. The noise spectra are 
precisely the measured quantity in experiment. In the 
above expression we have dchncd the spectral density of 
the random process y(t) to be, 



G v (uj) = lim — — 



T/2 



[y(t) - yje^dt 



-T/2 



(15) 



where y is the average of y(t) over time T. 

In summary, the explicit steps necessary to determine 
the stochastic dynamics of a cantilever array are: 

1. Perform the deterministic calculation. Given some 
arrangement of cantilevers in a fluid choose one 
cantilever and apply a step force some time in the 
distant past. This will eventually result in a fi- 
nite steady deflection of this cantilever, whereas the 
other cantilevers will have zero deflection. 

2. Remove the step force and calculate the deflections 
of the cantilevers as a function of time as they re- 
turn to equilibrium, for two cantilevers this is X\ (t) 
and X2(t) (see for example Fig.[|)]for the case of one 
cantilever). This deterministic calculation can be 
done analytically or numerically depending upon 
the complexity of the particular situation. 

3. Calculate the stochastic quantities. Using X\(t) 
and X2 (t) calculate the auto- and cross-correlations 
of the fluctuations from Eqs. I|llfl an d (|12|l . 
In essence, the cantilever deflections are merely 
rescaled to produce the stochastic response (this 
is illustrated by the two ordinate axes in Fig.^J). 

4. Use the auto- and cross-correlations to calculate the 
noise spectra given by Eqs. (jT3"|) and l(Ti|) . 



III. ANALYTICAL THEORY BASED ON AN 
OSCILLATING INFINITE CYLINDER 

The equipartition theorem permits the calculation of 
the mean square mode displacements of micron and 
nanoscale cantilevers in vacuum j^]. Sader |2^| further 
developed this approach to include the damping forces of 
a surrounding viscous fluid and solved the equations of 
elasticity in the thin beam limit coupled together with 



the Navier-Stokes equations. He used the approxima- 
tion that the effect of the fluid on each element of can- 
tilever is the same for an element of a long and slender 
cantilever moving with the same speed. In the limit of 
an infinitely long cylinder the flow field around the can- 
tilever can be assumed to be the flow field around an 
infinite two-dimensional cylinder to a good approxima- 
tion [2 St |29j. In addition, because the forces on each 
element of cantilever due to the fluid depend only on the 
velocity of that element, and not otherwise on its posi- 
tion in the cantilever, the mode structure of the damped 
cantilever is unchanged from the undamped limit. We 
follow this approach here and use the infinite cylinder 
approximation for the flow field. However we do not as- 
sume that the fluid noise is white as was done in |28|| but 
use the correct spectral density of the Brownian force 
given by Gj? (w) in Eq. i|37|) . The frequency dependence 
of the Brownian noise Gf (w) is determined by the fluid 
damping 7/(w) (illustrated in Fig. |21 and discussed fur- 
ther below). 

To connect with the previous work we note that the 
susceptibility x(t) connecting the tip displacement to an 
applied force is given by the cantilever response to a unit 
impulse of force, which is the derivative of the step force 
response given by Eq. (fTT|) so that 



(16) 



The Fourier transform of x(t) is the response to a deter- 
ministic sinusoidal force and the imaginary part — 
Im{x(w)}, where Im{} indicates the imaginary compo- 
nent is, 

f°° d 

X» = -0J j t (zi(O)zx(t)) sin(cjt)dt, (17) 

where we have adopted the Fourier transform convention 
given by, 



(18) 
(19) 



x(lu) = / x(t)e luJt dt 
J —00 

1 f°° 

x{t) = — / x(uj)e- lbJt du>, 



and i = y — 1. Integrating Eq. 1|17|) by parts and using 
the definition of the noise spectrum Eq. I|13|l yields, 



G u (u) 



x H- 



(20) 



The susceptibility x(t) can be determined by solving for 
the deterministic cantilever response to a force impulse. 

In the following we present results only for the funda- 
mental mode; however higher modes can be included if 
desired. The equation governing the cantilever motion is 
then. 



m e Xi + kXi = F f + Fi, 



(21) 
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where the amplitude of the mode displacement is char- 
acterized by the tip displacement Xi, m e is the effective 
mass of the cantilever in vacuum, Ff is the force acting 
on the cantilever due to the fluid, and Fi — 6(t) is the 
force impulse (8(t) is the Dirac delta). The effective mass 
of the cantilever m e is chosen to yield the same kinetic 
energy as in the cantilever mode, and is related to the 
cantilever mass by, 



am c , 



(22) 



where m c is the actual cantilever mass and for the funda- 
mental mode of oscillation for a beam a = 0.243. Taking 
the Fourier transform of this equation yields, 



-m P u> 



+ k)X l 



F 



f 



1. 



(23) 



The force from the fluid can be written in the form, 

Ff = m cy i >e u> 2 r(u))x, (24) 



where 



m cv i,e = otm C yi = apf 



f- 

V4 



w 2 L 



(25) 



Here m cy i is the effective mass of a fluid cylinder of ra- 
dius w/2 where pf is the fluid density, w is the cantilever 
width, and L is the cantilever length. Again the pref- 
actor a — 0.243 is to take into account the variation of 
the fluid force along the cantilever due to the varying 
velocity given by the mode structure. The fluid loading 
and damping are captured by the hydrodynamic function 
r(w) given by, 



r(w) = i + 



UK^-WiR) 



iRK Q (-iViR) 



(26) 



where K\ and Kq are Bessel functions jsfj. In this def- 
inition the frequency dependence on the right hand side 
appears through the frequency dependent Strouhal num- 
ber R. 

The cantilever is loaded by the fluid which can be char- 
acterized by an effective mass, m/, larger than m e that 
takes into account the fluid mass that is also being moved. 
The fluid also damps the motion of the cantilever which 
can be expressed as an effective damping jf. Relations 
for rrif and 7^ can be found by expanding T(ui) into its 
real and imaginary parts in Eq. (|23() and rearranging such 
that, 



to give, 



and, 



-mf(Lu)uj 2 Xi - iujj f (uj)Xi + kX t = 1 (27) 
mf = am c (l + T r') (28) 
7/ = am cy iLoT" , (29) 



where V and T" are the real and imaginary parts of T, 
respectively. T is the mass loading parameter, which is 
the ratio of the mass of a cylinder of fluid with radius 
w/2 to the actual mass of the cantilever, and is given by, 



T>cyl 



IT pfl 



4 p c h 



(30) 



where p c is the density of the cantilever and h is the 
cantilever thickness. 

It is evident from Eqs. and l|2*9"|l that both the fluid 
loaded mass of the cantilever and the fluidic damping are 
functions of frequency. The ratio of the mass of the fluid 
loaded cantilever to the effective mass of the cantilever 
in vacuum, m e , as a function of frequency is given by 
m//m e = 1+TqT' (R). The factor To is a constant for any 
particular choice of cantilever and fluid combination. The 
frequency dependence of the added mass is then given by 
T'(R) and is illustrated by the solid line in Fig.[3]and uses 
the left ordinate axis. The added mass increases rapidly 
as the frequency is decreased. Over the range of four or- 
ders of magnitude in frequency the added mass is seen to 
change by a factor of approximately 25. The frequency 
dependence of the fluid damping is given by loT" (uj) and 
is shown by the dashed line in Fig. [21 using the right or- 
dinate axis. There is a weaker frequency dependence for 
the fluid damping when compared to the mass loading. 
The fluid damping decreases as the frequency decreases 
and over 4 orders of magnitude in frequency the damping 
changes by a factor of approximately 7. 




FIG. 2: The frequency dependence of the fluid loaded mass 
and fluid damping for an oscillating cantilever in a viscous 
fluid, (solid line) T'(_R) illustrates the frequency dependence 
of the fluid loaded mass, (dashed line) RT"(R) illustrates the 
frequency dependence of the fluid damping. 

Solving for the cantilever response Xi(uj) and taking 
the imaginary part yields, 

UJJf 



X" 



(—rrifUj 2 + k 2 ) 2 — (^7/) 2 



(31) 
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Using this result in Eq. I|20|) and rearranging gives the 
desired result for the spectral density of the stochastic 
fluctuations in cantilever displacement, 



Ak B T 1 

i 

k cj 



(32) 



T u)Ti(R u)) 



(1 - £ 2 (1 + T T r (R u))) 2 + (^ToTiiRou)) 2 



where u> — lo/lo is a nondimcnsional reduced frequency, 
and i?o is the Strouhal number evaluated at the resonant 
frequency in vacuum given by, 



4 1] 



(33) 



We emphasize that Eq. I|32[l is not the same as in previous 
work [28| which failed to include the frequency dependent 
of the Brownian force. 

To understand the difference from the previous work 
we now connect Eq. I|10|) with the calculation in terms 
of a fluctuating force with spectral density Gf{u). The 
fluctuating displacement can be described in terms of the 
response to this force through the susceptibility 



GnM = IxMI'GfM. 

Inserting Eq. (|20|) into this expression yields, 



(34) 



(35) 



It is useful at this point to introduce the impedence 
Z(t) = F/v where v is the velocity of the cantilever tip 
(x\(t)), it is then straight forward to show, 



Z(co) 



1 



(36) 



Inserting this into Eq. 113511 and rearranging yields the 
desired expression of the fluctuation dissipation theorem, 



G F (u) = Ak B Tj f (u>) 



(37) 



where 7/(w) = Re |z(cj)| is the resistance or dissipation 

and Rc{} indicates the real part. For the case of oscillat- 
ing cantilevers in fluid 7/ is the effective fluid damping. 
It is evident from the frequency dependence of the damp- 
ing that the fluctuating force described by Eq. (|37Jl is not 
white noise as considered previously j2g (see Fig. [21 for 
the variation of 7/ with frequency). 

The spectral density of fluctuations in cantilever dis- 
placement G11 (a;) can also be determined by directly 
solving the governing stochastic equation of motion, 



m e x\ + kx\ = Ff + Fb, 



(38) 



where Fb is the random force due to Brownian motion. 
Taking the Fourier transform of this equation yields, 



(-m e uj 2 + k) x = F f + F E 



(39) 



Solving for the magnitude of the cantilever response and 
using Eq. (|37|l for the spectral density of the Brownian 
force again yields Eq. H32|) . Note that integrating this 
result automatically leads to the equipartition result, 



2tt 



Gii(oj)duj 



k ■ 



(40) 



We would like to point out that it is not possible to start 
with the equipartition result given by Eq. (|40|l and to 
then calculate the spectral properties of the fluctuations, 
one must first start with the fluctuation-dissipation result 
given by Eq. H37|) as done here. 

The noise spectrum of the cantilever in fluid may be 
used to calibrate an AFM. A simple way to do this is 
to extract the effective spring constant from the the fre- 
quency u> f giving the maximum of the noise intensity. In 
previous work |28l | the fluctuations of the cantilever have 
been estimated without the frequency dependence of the 
numerator in Eq. (|32|l . leading to some inaccuracy in the 
estimate of the spring constant. Once w/ is known, the 
quality factor Q can be estimated from, 



Q 



uorrif 
If 



j_ 

To 



T r (R) 



(41) 



where Eqs. (|28l) and (|29|l have been used for the mass 
and damping. This expression for Q is only approximate 
again because the frequency response of the cantilever is 
not precisely like that of simple harmonic oscillator, since 
both mf and 7/ are frequency dependent. 

The error incurred by neglecting the frequency depen- 
dence in the numerator in Eq. (l-il'li in fitting the spectrum 
of a Brownian driven cantilever is illustrated in Fig. as 
a function of Rq and To. In Fig. [3ju;| is the frequency 
at the maximum of the noise spectrum when the fre- 
quency dependence has not been included (specifically, 
in Eq. I|43|l Gi?(w) has been assumed constant) and ujf is 
the correct value of the frequency using the fluctuation- 
dissipation theorem as discussed in Section ITTT1 The re- 
sults indicate that as the Strouhal number increases the 
error in this assumption becomes smaller. In particu- 
lar, when the frequency dependence is neglected the ap- 
proximate theory always underpredicts ujf. This can be 
understood in light of Fig. where it is clear that the 
magnitude of the damping decreases as the frequency of 
oscillation is reduced. When this decrease is not included 
the predicted value of ujf will be unnecessarily reduced. 
Although the fluid induced damping slowly decreases as 
the oscillation frequency is reduced it is important to 
note that the fluid loaded mass increases rapidly and it 
is this interplay which results in the small Q associated 
with oscillating micron and nanoscale cantilevers. 

Figure 01 summarizes the theoretical results and allows 
for easy determination of the stochastic dynamics of a 
single cantilever of arbitrary geometry placed in an ar- 
bitrary viscous fluid. Given a particular cantilever and 
fluid combination the procedure is: 
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T 



FIG. 3: The error in predicting the frequency at the max- 
imum of the noise spectrum u)f/uf. oj} is the predicted 
frequency at the maximum of the noise spectrum when the 
frequency dependence of the Brownian force is neglected as 
done in previous work |2SJ an d 10 f ls the correct value given 
by the results of Section lllll Results are presented as a 
function of mass loading To and the Strouhal number Ro 
based upon the cantilever's resonant frequency in vacuum. 
The curves are for R = 128, 64, 32, 16, 4, 1, 0.25, 0.0625 where 
Ro = 0.0625, 1, 128 are labelled and the remaining curves are 
in sequential order. 

1. Determine the cantilever spring constant k, reso- 
nant frequency in vacuum, ujq, and characteristic 
length scale. These can be experimental measure- 
ments or theoretical calculations. The character- 
istic length is a measure of the characteristic half- 
width of the cantilever. 

2. Determine Ro from Eq. (J2SI and T from Eq. (J30). 

3. Use Ro and To to determine ujf/u) from Fig. Ufa) 
and Q from Fig.^b). 

In the case that the cantilever is driven externally the 
response can be found in a similar manner. If the can- 
tilever is driven by a force Fd(t) the equation of motion 
becomes, 

m e x + kx = F f + F d . (42) 

If the driving force is Fd — Ao sin(wdt) where Aq is a 
constant and uid is the driving frequency the amplitude 
of the response is given by 

|iM| 2 =(^f) 2 . (43) 

1 

(1 - u\ (1 + T T r (R u d ))) 2 + (ujToT^RoUd)) 2 ] ' 
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FIG. 4: Panel (a) the reduced frequency u> of oscillation for a 
cantilever undergoing stochastic oscillations in a viscous fluid 
as a function of mass loading To and the Strouhal number Ro 
based upon the cantilever's resonant frequency in vacuum. 
Panel (b), the quality factor Q as a function of Ro and Tb 
as given by Eq. 141H . In each panel the largest value of Ro 
is labelled and each successive curve represents Ro/2, i.e. in 
panel (a) the last curve is for Ro = 0.0625 and in panel (b) 
the last curve is for Ro = 1/2. 



where uj d — ^d/^a is the reduced driving frequency. We 
emphasize that Eqs. H42|) and H43fl ar< 3 for externally 
driven cantilevers and neglect Brownian noise. Eqs. I)43|) 
and (|32|l are similar in that they share a common denom- 
inator, whereas the additional frequency dependence in 
the numerator of Eq. I|32|l is from the frequency depen- 
dence of the Brownian noise. 



IV. MICRON SCALE BEAM IN FLUID 
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Q m f /m e Uf/wo 7/ (kg/s) R f 



The stochastic dynamics of a micron scale cantilever, 
similar to what is commonly used in atomic force mi- 
croscopy, are now quantified (se e Fig. [SJ. The beam is 
chosen from Chon and Sader [31| where both theoretical 
and experimental values are present for comparison. The 
beam properties are summarized in Table [H using clas- 



wall 




FIG. 5: Schematic of a simple cantilevered beam of length L, 
width w, and height h. 



fo 



Ro T 



197/im 29^m 2pm 71.56KHz 109.7 4.89 

TABLE I: The beam geometry: length L, width w, thickness 
h, resonant frequency in vacuum fo, Strouhal number based 
on the vacuum resonant frequency Ro and the mass loading 
factor T . The beam is silicon with E = 1.74 x 10 n N/m 2 , 
p c = 2320Kg/m 3 and the fluid is water with (p f = 997Kg/m 3 , 
77 = 8.59 x 10 _4 Kg/ms). In the deterministic numerical simu- 
lations the magnitude of the step force applied is To = 26nN. 



sical beam theory [3^ 
predicted to be, 



the equivalent spring constant k is 



3EI 



(44) 



which yields a value of k = 1.3N/m and the resonant 
frequency of the cantilever in vacuum is 



fo 




(45) 



which yields f = 71.56KHz where C\ = 1.8751. Using 
Eq. \.YM) the Strouhal number in vacuum is Ro = 109.7 
and using Eq. (|3Tj|l the mass loading factor is T — 4.89. 

We now use the thermodynamic approach to calcu- 
late the stochastic dynamics of the cantilever. As pre- 
viously stated this can be accomplished in a straightfor- 
ward manner by determining the deterministic response 



(1) 

(2) 



3.24 
2.93 



8.16 
8.05 



0.34 
0.35 



5.07 x 10" 7 
6.91 x 10" 7 



37.3 
38.7 



TABLE II: The stochastic dynamics of a single micron scale 
cantilever in fluid. (1) Results based on the analytical predic- 
tions of Section UTTl The quality factor Q is calculated from 
the approximation given by Eq. (gTJ. (2) Results from fi- 
nite element numerical simulations using the thermodynamic 
approach of Section [HI The numerical results are fit to the 
simple harmonic oscillator response given by Eq. 14711 to de- 
termine the listed diagnostics. 



of the cantilever to the removal of a step force applied to 
the tip. We do this using finite element numerical sim- 
ulations of the full three dimensional, time-dependent, 
fluid-solid interaction problem (algorithm discussed else- 
where H3,I3)- The simulation is initiated with the re- 
moval of a step force applied to the tip of the cantilever 
and the deterministic dynamics of the beam Xi(t) are 
shown by the solid line in Fig. fusing the right ordinate 
axis. 

In order to calculate commonly used diagnostics, such 
as ujf and Q, the deterministic cantilever deflection X\(t) 
was fit to the deflections of an underdamped simple har- 
monic oscillator (valid for Q > 1/2) given by the equation 
of motion, 

mfX^+jfX^ + kXiit) =0, (46) 
where X^O) = Fi/k and Xi(0) = 0. The solution is, 
Xl (t) = ^ e -/*/2Q ^ cos(a/t) + sin(c / t) ^ (47) 

where, 



1 



1 



(48) 



and Wf = \Jkjmf. Using a nonlinear least squares curve 
fit algorithm [35| the numerical results are fit to the 
model to yield the values of Q, u>f, 7/, and m/ shown in 
Table ITU The curve fit is nearly indistinguishable from 
the numerical simulation results for Xi{t) and is shown 
by the dashed line in Fig. EjJ 

Using Xi(t) in Eq. (|llfl yields the autocorrelation 
of the equilibrium deflections, (xi(O)xi(t)) , which are 
shown in Fig. El using the left ordinate axis. The noise 
spectrum is given by Eq. (|13fl and is shown in compari- 
son with predictions of the infinite-beam approximation 
given by Eq. (|32|l in Fig. The numerical simulations 
contain all of the oscillation modes as indicated by the 
presence of the second mode in the figure. In the infinite 
beam approximation the different modes behave indepen- 
dently and the calculation could be extended to include 
higher modes if desired. In Table [H] the analytical re- 
sults of Section IIIII are compared with results using the 
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full thermodynamic approach and deterministic finite el- 
ement numerical simulations. In general, from the table 
and the figures, we see that the differences between the 
simple model and the full calculations are small for this 
shape cantilever. The experimentally measured value of 
the resonant frequency in water is LUf/ujo = 0.36 [3ll | 
which agrees well with both the analytical results and 
the numerical predictions using the thermodynamic ap- 
proach. However, Table [H] suggests that the analytical 
results based upon the oscillating infinite cylinder model 
under predict the amount of damping jf. The noise spec- 
tra show differences that are small, but may be significant 
for quantitative applications such as calibration. 

In all of the deterministic finite clement numerical sim- 
ulations performed we were careful to ensure that the 
bounding no-slip surfaces of the computational domain 
did not affect the results. Comparisons with numerical 
simulations performed with larger domains did not sig- 
nificantly alter the results. 



0.003 




t (ms) 



FIG. 6: The dynamics of a micron scale cantilever immersed 
in water, (solid line) Deterministic finite element numeri- 
cal simulations using the thermodynamic approach, (dashed 
line) A simple harmonic oscillator curve fit to the data. The 
left ordinate yields the autocorrelations of the fluctuations in 
cantilever displacement and the right ordinate yields the de- 
terministic cantilever dynamics in response to the removal of 
a step force. 



V. AN ARRAY OF NANOSCALE 
CANTILEVERS 

As the cantilever dimensions become smaller the effec- 
tive cantilever spring constant decreases while the res- 
onant frequency increases. This favorable combination 
potentially provides access to the biologically important 
parameter regime characterized by 10's of piconewtons 
with microsecond scale time resolution, a range that is 



t 1 1 1 1 1 1 1 1 1 1 1 1 r 




FIG. 7: The noise spectrum Gu(u) for a micron scale can- 
tilever in water, (solid line) Deterministic finite element 
numerical simulations using the thermodynamic approach, 
(dashed line) Approximate analytical theory for the funda- 
mental mode only. 



difficult to reach using other methods. This has led to 
the development of nanoscale cantilevers [lil fl8L f36j | . 

In what follows we quantify the stochastic dynamics of 
two adjacent nanoscale cantilevers immersed in water. A 
schematic of the nanoscale cantilever under consideration 
here is shown in Fig. El 16lll8|. This is an experimentally 
motivated cantilever whose paddle shaped geometry de- 
creases the effective spring constant while localizing the 
strain for piezoresistive measurement. The configuration 
of the cantilever array is shown in Fig. ^ with two oppos- 
ing cantilevers separated by a distance s. This configura- 
tion was chosen because of its experimental accessibility 
as well as its potential use in single molecule measure- 
ments with the tethering of a target biomolecule between 
the two cantilevers. In the work presented here we build 
a baseline understanding of the cantilever dynamics in 
the absence of target biomolecules. Further considera- 
tion of the dynamics caused by a tethered biomolecule 
is beyond the scope of the present paper. The dimen- 
sions of the nanoscale cantilevers investigated here are 
summarized in Table ITTll 

As discussed previously, the physical properties de- 
scribing the dynamics of the cantilever in vacuum are 
important parameters in determining the stochastic dy- 
namics of the cantilevers. These properties are summa- 
rized in Table ITVl where the values have been determined 
from finite element simulations of the elastic cantilever 
structure in the absence of the surrounding fluid. It is 
important to emphasize that although these calculations 
can sometimes be performed analytically, this is not the 
case for many geometries of interest. However, even for 
cantilevers with very complex geometries, the vacuum 
based finite element calculations are straight forward. 
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L w h L\ b 



3pm lOOnm 30nm 0.6pm 33nm 

TABLE III: Geometry of the nanoscale cantilever (see Fig.[SJ. 
The nanoscale cantilever is composed of silicon (E = 1.25 x 
10 n N/m 2 , p c = 2330kg/m 3 ) and the fluid is water (p f = 
997kg/m 3 , r) = 8.67 x 10~ 4 kg/ms). The magnitude of the 
step force applied to the tip of the cantilever is Fo = 0.75pN. 

k coo Ro To 



8.7 mN/m 37.46 x 10 6 rads/s 0.11 1.12 

TABLE IV: The cantilever spring constant k, resonant fre- 
quency in vacuum ljo , the Strouhal number R$ , and the mass 
loading factor To. ljq and k were determined from numerical 
simulations of the nanoscale cantilever in vacuum. 



The stochastic dynamics of the cantilever array are de- 
termined using the thermodynamic approach discussed 
in Section [n] with deterministic finite element numeri- 
cal simulations. A series of numerical simulations have 
been performed to determine the variation in the can- 
tilever dynamics as a function of the separation distance 
between the two cantilevers in the array. In particular we 
have explored s/h = 1, 2, 3, 4, 5 where s is the cantilever 
separation and h = 30nm is the cantilever thickness. The 
numerical simulations are initiated by the removal of a 
step force applied to the tip of one of the cantilevers (in 
Fig. n this is the cantilever on the left) while the other 
cantilever is initially undeflected and at mechanical equi- 
librium. 

The deterministic dynamics of the nanoscale cantilever 
for which the step force has been removed, Xi(t), is 
shown by the solid line in Fig. [5] using the right ordi- 
nate axis. The deflection -Xi(t) was not significantly af- 
fected by the presence of the second cantilever for the 
separations considered here. As with the micron scale 



Wall 




w z 



FIG. 8: Schematic of an experimentally motivated nanoscale 
cantilever geometry of length L, width w, and height h. 



cantilever, we again fit the result for X\ (t) from the finite 
element numerical simulation to the solution of the sim- 
ple harmonic oscillator equation given by Eq. (|46|l . For 
the nanoscale cantilever the dynamics are overdamped 
(Q < 1/2) and the deterministic deflection predicted by 
the simple harmonic oscillator approximation is, 

where, 

^=^(-^ ± ^5^)- (50) 

The agreement with the curve fit is very good and uif 
and Q are shown in Table [YJ Also shown in Table El 
are the theoretical predictions based upon the infinite 
cylinder approximation. For the analytical calculations 
of the damping "/f from Eq. I|29|l the fundamental mode 
of the nanoscale cantilever has been modelled as that of 
a hinge where all of the bending occurs in the short legs 
near the base. This is guided by numerical calculations 
showing the strain localized in the short legs. In this case 
we find the mass factor is a — 1/3. (Note that a is not 
needed to determine the other quantities in Table IVl) . It 
is clear from Table IVl that the infinite cylinder approxi- 
mation is no longer valid. The analytical prediction for 
the added mass is an order of magnitude too large. This 
suggests that three-dimensional flow effects become sig- 
nificant for the shape of cantilever under consideration 
here. For example there would be flow around the tip 
of the cantilever, as well as the flow through the open 
region near the base of the cantilever. In addition, the 
predicted frequency shift is an order of magnitude too 
small. 

Q m f /m e LUf/ujQ 7/ (kg/s) R f 



(1) 0.265 106.3 0.032 2.41 x 10~ 9 0.0036 

(2) 0.327 18.16 0.236 1.01 x 10~ 9 0.026 

TABLE V: The stochastic dynamics of a nanoscale cantilever 
in fluid. (1) Results based on the analytical predictions of Sec- 
tion ^H] The quality Q is calculated from the approximation 
given by Eq. (I4H . (2) Results from finite element numerical 
simulations using the thermodynamic approach of Section ITT1 
The numerical results are fit to the simple harmonic oscillator 
response given by Eq. 1491 to determine the listed diagnostics. 

Inserting Xi(t) into Eq. (|llfl yields the autocorrela- 
tion of the equilibrium fluctuations in cantilever displace- 
ment, (xi(0)xi(t)} . These are also shown by the solid 
line in Fig. [U] using the ordinate axis on the left. The 
noise spectrum, Gu(w), is found by inserting the auto- 
correlation in Eq. (|13fl and the result is shown by the 
solid line in Fig. Fig- EH a l S0 presents a comparison 
of the noise spectrum Gn(w) given from the theoretical 
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FIG. 9: Comparison between the model theory and the 
full numerical simulation for the autocorrelation of equilib- 
rium fluctuations (right ordinate axis) and the determinis- 
tic cantilever deflection relative to the initial deflection of 
Xq — Xi (t — 0) (left ordinate axis) for a single nanoscale can- 
tilever in fluid: (solid line) results from finite element numer- 
ical simulations using the thermodynamic approach; (dashed 
line) results based upon the approximate analytical theory of 
Section ITTT1 



predictions of Eq. (|32[l and the results from the numeri- 
cal simulations using the thermodynamic approach. The 
noise spectra are different for small frequencies, in partic- 
ular the model result has a peak at finite frequency not 
seen in the full calculation. The inverse cosine transform 
of the noise spectrum from the infinite-cylinder model 
yields the deflection of the cantilever as a function of 
time shown by the dashed line in Fig. [5] using the right 
ordinate axis. The model prediction for the cantilever 
deflection exhibits negative values that are not seen in 
the simulations suggesting that the response of the over- 
damped cantilever is not precisely that of a simple har- 
monic oscillator. 

We now turn to the deterministic deflection of the sec- 
ond cantilever, X2(t), after removing the force on the first 
cantilever, evaluated using the full fluid-elasticity simu- 
lations. This is shown in Fig. ^Jusing the right ordinate 
axis. For close separations s/h < 1 the second cantilever 
exhibits negative deflections for all time. However, as 
the cantilever separation increases the second cantilever 
initially exhibits a positive deflection followed by a nega- 
tive deflection. The flow field around an oscillating object 
will have both potential and non-potential components. 
For an incompressible fluid the potential component is 
instantaneous whereas the non-potential component dif- 
fuses with a diffusion coefficient given by the kinematic 
viscosity v. Although the deterministic flow field around 
simple oscillating objects is well known the fluid coupled 
motion of multiple elastic objects is not well understood. 
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FIG. 10: The noise spectrum Gri(w) for a single nanoscale 
cantilever immersed in water, (solid line) Results from deter- 
ministic finite element simulations using the thermodynamic 
approach described in Section [H] (dashed line) Results from 
approximate analytical theory described in Section ITTT1 

The fluid dynamics resulting from the motion of two ad- 
jacent cantilevers will be discussed in a forthcoming ar- 
ticle H3- 

Using Eq. (|12[) the stochastic correlations between the 
two cantilevers can be determined from A^i). The cross- 
correlations of the equilibrium displacement fluctuations 
(xi(0)x2(t)} are shown in Fig. using the ordinate axis 
on the left. For small cantilever separations s/h < 1 the 
cross-correlations are anticorrelated for all time. This 
is in agreement with experimental measurements of the 
fluid correlations of two closely spaced micron scale beads 
in water [l^ . l38| . However, as the cantilever separation is 
increased the cross-correlations changes: for short times 
the cross-correlations are positively correlated whereas 
for larger times they are anti-correlated. 

The noise spectra Gi2(w) are found from (xi(0)x2(t)) 
by using Eq. Q14|). Figure IT51 shows the variation in the 
noise spectra as a function of cantilever separation. The 
noise spectra contain both positive and negative values 
and for each cantilever separation there is a frequency at 
which correlated noise vanishes, Gi2(w) = 0. This null 
point could be exploited by a measurement scheme to 
minimize the correlated noise. 

From these results it is possible to characterize the 
force sensitivity and time resolution of a correlation 
measurement technique using an array of closely spaced 
nanoscale cantilevers. An estimate of the force sensitiv- 
ity can be found using the auto and cross correlation 
functions as, 

F = fcllMO)^))^, (51) 
F s = fc|Kzi(0)z 2 (i))||V a 2 x . (52) 

Fq represents the approximate magnitude of the stochas- 
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FIG. 11: The cross-correlations in equilibrium fluctuations of 
cantilever deflections for two adjacent nanoscale cantilevers 
separated by distances s/h — 1, 2, 3, 4, 5 for h = 30nm found 
using deterministic finite element numerical simulations with 
the thermodynamic approach discussed in Section [HI The 
curve for s/h — 1 is labelled and the others follow in sequential 
order. The right ordinate axis illustrates the deterministic 
cantilever deflection as a function of time X2(t) scaled by the 
deflection of the adjacent cantilever, Xo = X\(t = 0). 



tic Brownian force acting on a single cantilever and F s 
the force induced by fluid correlations between two can- 
tilevers. For the case of a single nanoscale cantilever the 
maximum value occurs at t = and, using the equiparti- 
tion theorem this yields a value of F — y/ksTk = 6.0pN. 

In a cross-correlation measurement between two can- 
tilevers the Brownian noise felt by the two individual 
cantilevers is uncorrelated and does not contribute. This 
leaves only the correlations due to the viscous fluidic cou- 
pling and F s represents the approximate magnitude of 
the force due to this hydrodynamic coupling. For exam- 
ple, from Fig. ^] for the case with the closest separa- 
tion s/h — 1, the magnitude of the maximum value of 
the cross-correlation is (xx(0)x2(t)) — 1.19 x 10 _2 nm 2 
which yields a force sensitivity F s = 0.95pN. Therefore 
F 8 /Fq w 1/6 at s/h — 1, indicating a 6-fold reduction in 
thermal noise over a single cantilever measurement. 

The variation in F s as a function of cantilever sepa- 
ration is shown in Fig. I13f ah The data is fit with an 
exponential curve, given by F s = 1.12 exp(— 0.18s/h) (in 
units of pN), shown in Fig. I13t a) as the solid line. For 
very small cantilever separations s/h < 1 the noise re- 
duction possible through the use of a correlated measure- 
ment technique is approximately 6-fold. As expected, the 
correlated noise decreases as the cantilever separation is 
increased: for a separation of s/h = 5 there is a 12-fold 
noise reduction. The ordinate axis on the right hand 
side of Fig. I13f a) shows F s /Fq to show this performance 
improvement directly. 

The characteristic time scales of the fluid correlated 
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FIG. 12: The noise spectrum Gi2(u>) for two nanoscale can- 
tilevers separated by distance s/h = 1, 2, 3, 4, 5 for h = 30nm 
found using deterministic finite element numerical simulations 
with the thermodynamic approach discussed in Section [H] 
The curve for s/h = 1 is labelled and the others follow in 
sequential order. 



motion as a function of cantilever separation are illus- 
trated in Fig. ICT b). In Fig. ICT b) circles represent the 
time at which the maximum value in the cross correla- 
tion occur and squares represent the time at which the 
fluid induced correlations vanish (i.e., the zero-crossing 
in Fig. Ill|) . As expected from the finite rate at which 
momentum diffuses away from an oscillating cantilever, 
both of these time scales increase with cantilever separa- 
tion. Using the ordinate axis on the right in Fig. 113T b) 
the ratio of these time scales to the period of a single os- 
cillation for the cantilever in vacuum to = Jq 1 are given. 



VI. CONCLUSIONS 

The stochastic dynamics of micron and nanoscale can- 
tilevers can be quantified for the precise conditions of ex- 
periment using straightforward deterministic calculations 
coupled with the fluctuation-dissipation theorem. This 
can be done for complex cantilever geometries includ- 
ing wall effects, and the stochastic correlated behavior 
of closely spaced cantilevers for which theoretical predic- 
tions currently are not available. Cantilever geometries 
that are long and thin are well described by approxi- 
mate analytical theory based upon the oscillation of an 
infinite cylinder in fluid. We have used the fluctuation- 
dissipation method to correct the previous results of this 
theory. For the short and wide nanocantilevers currently 
being proposed the infinite-cylinder approximation is no 
longer appropriate. The thermodynamic approach us- 
ing the fluctuation-dissipation theorem provides an im- 
portant means of developing a physical understanding 
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FIG. 13: The force sensitivity and time scales characterizing 
two adjacent nanoscale cantilevers found using deterministic 
finite element numerical simulations with the thermodynamic 
approach discussed in Section [IT] Panel (a), the maximum 
magnitude of the force induced by fluid-coupled correlations 
in cantilever displacement F s . The data is empirically fit with 
an exponential given by F s — 1.12 exp(— Q.18s/h) (in units of 
pN). The magnitude of the force felt by a single cantilever due 
to Brownian noise, Fo = 6.0 pN, is nearly 6 times larger than 
the correlated noise at close separations (s/h < 1). The right 
ordinate axis scales F s by Fo. Panel (b), characteristic time 
scales, the circles represent the time at which the maximum 
value of the cross correlation occurs, the squares represent 
the time at which the fluid induced correlations vanish. The 
ordinate axis on the right hand side illustrates the time scales 
when scaled by the period of oscillation of the cantilever in 
vacuum to = (cvo/2w)~ 1 . The solid lines represent a linear 
curve fit to the data. 
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